Simulating Quantum Magnetism with 
Correlated Non-Neutral Ion Plasmas 

John J. Bollinger, Joseph W. Britton and Brian C. Sawyer 

Time and Frequency Div., NIST, Boulder, CO 80305, USA 

Abstract. By employing forces that depend on the internal electronic state (or spin) of an atomic 
ion, the Coulomb potential energy of a strongly coupled array of ions can be modified in a 
spin-dependent way to mimic effective quantum spin Hamiltonians. Both ferromagnetic and anti- 
ferromagnetic interactions can be implemented. We use simple models to explain how the effective 
spin interactions are engineered with trapped-ion crystals. We summarize the type of effective spin 
interactions that can be readily generated, and discuss an experimental implementation using single- 
plane ion crystals in a Penning trap. 
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INTRODUCTION 

Currently, many important problems in physics, for example solid-state physics, are 
poorly understood because the underlying quantum mechanics is too complex for mean- 
ingful modeling. Frequently cited examples include correlated magnetic systems, such 
as spin liquids [1], and high-temperature (high-Tc) superconductivity, where we are 
not able to calculate the phase diagram or even know what phases exist in the Fermi- 
Hubbard model [2], a leading candidate for explaining high-Tc superconductivity. The 
source of this problem is the exponential increase in the difficulty of quantum many- 
body calculations as the number of bodies increases. Considering simple two-level quan- 
tum systems or spins, the complexity of the calculation doubles with each additional 
spin. This exponential increase in complexity with the number of spins makes a direct 
calculation of a general system of more than ~ 30 interacting spins intractable on current 
computers [3]. 

In an attempt to tackle this difficulty, a number of groups are following a suggestion 
put forth by Feynman [4] that it might be possible to engineer interactions between well 
controlled quantum components in order to mimic a quantum many-body system that 
is not understood. Through measurement of the well controlled quantum components, 
it may then be possible to obtain the solution, or at least acquire some information 
about the behavior of the currently intractable quantum many-body system. Current 
platforms (i.e., well controlled quantum systems) being used for quantum simulation 
include neutral atoms in optical lattices [5], superconducting circuits [6], and trapped 
ions [7-11]. Each quantum platform has particular strengths for engineering different 
types of many-body Hamiltonians. Neutral atoms in optical lattices are well suited 
for simulating many-body Hamiltonians where quantum statistics (Fermi or Bose) are 
important. Crystalline arrays of trapped ions are a promising platform for emulating 



quantum magnetic systems. Superconducting circuits provide exceptional flexibility in 
the design of interactions between elementary quantum building blocks. 

In the next section we use a simple model to illustrate how crystalline arrays of trapped 
ions can be used to emulate an interacting system of spin- 1/2 particles. Crystalline 
arrays of trapped ions are strongly coupled; the potential energy due to the Coulomb 
interaction is much larger than the ion thermal energy. The basic idea behind engineering 
quantum magnetic interactions is to generate an internal- state (or spin) dependence of 
the Coulomb interaction energy, through the application of forces that depend on an ion's 
internal state. In the final section we describe our experimental set-up for implementing 
quantum magnetic interactions on single -plane triangular arrays of ions stored in a 
Penning trap [11]. 

QUANTUM MAGNETIC INTERACTIONS THROUGH 
SPIN-DEPENDENT FORCES 

Ions have a variety of different internal states. For example, the ground and excited 
electronic configurations of an ion give rise to different internal states. Each of these 
configurations is typically split into a number of different hyperfine states corresponding 
to different orientations of the electronic spin and nuclear moment. With well developed 
techniques from atomic physics, it is frequently possible to isolate and control two 
different internal atomic levels — a two-level system or qubit — with high fidelity. This 
two-level system acts as an effective spin- 1/2 [12] and we denote these two levels with 
the notation | |), | !}. In our experimental set-up (described later), the effective spin- 
1/2 is the valence electron spin-flip transition in the ground state of Be + in the 4.46 T 
magnetic field of the Penning trap. In this case the isolated two-level system is very close 
to a real spin- 1/2 particle, but this does not need to be the case. Because the distance d 
between ions in a trapped ion crystal is large (typically d ~many micrometers), the 
interaction energy between two neighboring ions in the crystal lattice is simply the 
spin-independent Coulomb interaction e 2 / (4ne d). For the valence-electron spin qubit 
discussed here, the magnetic dipole interaction energy at a separation of d ~ 10 fim is 
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f|^f ~ 8.6 x 10~ 39 J ~ h(12 /iHz), which is too weak to be of any consequence. 

Larger spin-spin interactions with trapped ions are obtained through the application of 
spin-dependent forces. A detailed description can get complicated, but the basic idea can 
be illustrated for the simple case of a two-ion crystal, as shown in Fig. 1. Suppose the 
ions are confined by a harmonic trapping potential characterized by a frequency CO z in 
the direction transverse to their separation . With two ions there are 4 different spin states 
shown in Fig. 1(a) that need to be considered: two ferromagnetic states (| t)| t)> 1 1)1 1)) 
and two antiferromagnetic states (| t)| 4-)> I 4-)| t))- Now suppose we adiabatically turn 
on a spin-dependent force in a direction transverse to the separation of the two spins (the 
z direction in Fig. 1). The spin-dependent force could be, for example, due to an applied 
inhomogeneous magnetic field or an optical dipole force. For simplicity we assume that 
the force on |f) is equal in magnitude but opposite in sign to the force F\ on |!}. If, for 
the moment, we neglect the Coulomb repulsion between the ions, then each spin |t) (ID) 
state is displaced from the trap center by a distance Az = Fo/(m(0 2 ) (Az = — Fq/ (m(D 2 )) . 
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FIGURE 1. (a) Simple example of two trapped ions separated by a distance d. In a direction normal 
to their separation the ions reside in a harmonic potential characterized by frequency co z . The four 
different spin states are shown, (b) Application of a static, spin-dependent force produces spin-dependent 
displacements of the ions. If the magnitude |Az| of the displacements are identical for all the spin states, 
the change in the Coulomb potential energy looks like an effective antiferromagnetic interaction. 



Because all ions are displaced from the trap center by the same distance, the increase in 
the external trapping potential energy bJij rap is the same for all spin states. However, 
for the anti-ferromagnetic states (|t)|4) 5 |4) |t))> th e Coulomb interaction energy is 
reduced because of the increased distance between the ions. This change in the Coulomb 
interaction energy mimics an antiferromagnetic interaction between the two ion qubits. 
This illustrates, at its most basic level, the idea behind generating effective spin-spin 
interactions (and quantum gates [13, 14]) with trapped ions. 

However, an important detail has been left out of Fig. 1. In particular, if we include 
the Coulomb repulsion between the ions, then the displacement of an ion from the trap 
center depends on the spin state of the other ion. For the antiferromagnetic states, the 




FIGURE 2. For the antiferromagnetic states, the Coulomb repulsion between the ions produces a larger 
displacement of the ions from trap center and an increase in the external trapping potential energy relative 
to the ferromagnetic states. The total change in energy mimics an ferromagnetic interaction. 



Coulomb repulsion enhances the effect of the spin-dependent force, producing a larger 
displacement of the ions from trap center and an increase in the external trapping po- 
tential energy relative to the ferromagnetic states. As sketched in Fig. 2, for small dis- 
placements (Az «d) the increase in the trapping potential energy is twice the decrease 
in the Coulomb interaction energy, with the net result that a static, spin-dependent force 
actually produces an effective ferromagnetic Ising interaction of strength 
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An engineered ferromagnetic interaction can be used, along with a transverse field, to 
study a paramagnetic-to-ferromagnetic quantum phase transition [8]. However, for the 
triangular lattice discussed in the next section, an antiferromagnetic interaction would 
be particulary interesting, as it enables studies of magnetic frustration and possible spin- 
liquid behavior. 

An antiferromagnetic interaction can be engineered through application of an oscillat- 
ing spin-dependent force F^(t) = —F\(t) = Fq sin /i^f . Once again, consider slowly ramp- 
ing up the amplitude of the spin-dependent force. We are interested in the steady-state 
modification to the energy of the system due to the oscillating spin-dependent force. If 
the frequency of the applied spin-dependent force is greater than the trapping frequency 



()1r > C0 Z ), the driven oscillation of the ions is 180° out of phase with the applied force, 
the exact opposite of that discussed above for the \1r = 0. For the antiferromagnetic 
states the Coulomb repulsion between the ions now opposes the spin-dependent force, 
resulting in a decrease in the driven amplitude of the antiferromagnetic states (relative 
to the ferromagnetic states). For \x r > C0 Z the energy of the driven motion is therefore 
greater for the ferromagnetic states than for the antiferromagnetic states, resulting in an 
antiferromagnetic Ising interaction. 

These simple examples illustrate the basic idea behind engineering spin-spin inter- 
actions with trapped ion crystals. Through the application of spin-dependent forces, 
the total energy of the crystal can be modified in a way that depends on the internal 
"spin-state" of the ions, generating an effective spin-spin interaction. The sign of the 
interaction can be tuned by employing a sinusoidally varying spin-dependent force: for 
Hr > C0 Z the effective interaction is antiferromagnetic; for }i r «(o z the interaction is fer- 
romagnetic. More complicated time-dependent forces can generate arbitrary spin-spin 
interactions [15]. 

ENGINEERED ISING INTERACTIONS ON SINGLE-PLANE 

CRYSTALS 

In this section we discuss our experimental set-up for generating Ising interactions 
between several hundred 9 Be + ions localized in single-plane triangular arrays con- 
fined by a Penning trap. Axial confinement in the direction of the magnetic field (£- 
direction, Bq=AA6 T) is due to the electric field from a quadrupolar trap potential 
(j)r = jniB e CO z (z 2 — p 2 /2), where typically co z ~ 2^x795 kHz. Rotation at frequency 
C0 r (about z) produces a radial restoring potential due to the velocity-dependent Lorentz 
force. co r is precisely controlled with an external rotating quadrupole potential (a "rotat- 
ing wall") [16]. For a weak rotating wall, the trapping potential in the rotating frame is 
q${p,z) ^ ^m Be (O z (z 2 + /3p 2 ), where j8 = (O r O)- 2 (Q c - (O r ) - \. For 100 < N < 350, 
we set (O r ~ 2k x 45 kHz so that the radial confinement is weak enough that a cloud of 
ions relaxes into a single 2D plane (/3 << 1). When the ions' motional degrees of free- 
dom are Doppler laser cooled (< 1 mK), the ions naturally form a 2D Coulomb crystal on 
a triangular lattice, which is the geometry that minimizes their mutual Coulomb poten- 
tial energy [17]. See Fig. 3(a) for an image of the triangular lattice taken in the rotating 
frame of the ions. 

Reference [18] gives a detailed discussion of our spin initialization, control, and mea- 
surement capabilities with planar ion arrays in Penning traps. Here we briefly summarize 
some of that discussion. Figure 3(b) shows the relevant 9 Be + energy levels. We use the 
valence-electron spin states parallel |t) = \mj = + \) and antiparallel ||) = \mj = — \) 
to the applied magnetic field of the Penning trap as the two-level system or spin. In 
the 4.46 T magnetic field of the trap, these levels are split by approximately 124 GHz. 
Spins in the |!) state are efficiently optically pumped to the \\) state by a laser tuned 
to the \\) — > \2P 3 / 2 , m J = +2) transition. The repump beam and the main Doppler laser 
cooling beam are directed along the magnetic field (z-axis). In addition, a weak Doppler 
laser cooling beam (40 /im waist) directed perpendicularly to the z-axis directly Doppler 
cools the perpendicular degrees of freedom. A typical experimental cycle starts with 10 




FIGURE 3. (a) Top-view resonance fluorescence image showing the center region of an ion crystal 
captured in the ion's rest frame [19]. The nearest neighbor ion spacing is ^20 jim. (b) Relevant energy 
levels of 9 Be + at Bq=4A6 T (not drawn to scale). The 9 Be + nucleus has spin 7 = 3/2. We show only 
mj = 3/2 levels that are prepared experimentally through optical pumping. The 2 5 , 1 /2 — 2 P3/2 transition 
wavelength is ~3 13 nm. The qubit splitting AEq ~ 124 GHz. A low-phase-noise microwave source at 124 
GHz provides full global control over spins, (c) The optical dipole force (ODF) interaction is due to a 
pair of beams (derived from the same laser) with relative detuning )1r. The beams cross with an angle of 
Or ~ 5° generating an optical dipole force with an effective wavelength of Xr = 2n/\Ak\ ps 3.7 jim. 



ms to 20 ms of combined Doppler laser cooling and repumping. The fidelity of the |t) 
state preparation is estimated to be very high (99.9 %). 

Low-phase-noise microwave radiation at 124 GHz is used to rotate the spins. The 
microwave source consists of an agile 15.5 GHz source followed by an amplifier mul- 
tiplier chain with 150 mW output power at 124 GHz. The microwaves are transported 
to the ions down the bore of the magnet with a rigid waveguide and directed onto the 
ions with a horn located between the ring and endcap electrodes of the trap. We obtain 
^-pulses (the duration of time required to coherently drive |t) to of 70 /is duration. 
The measured spin-echo coherence duration (T2) is ~100 ms. 

At the end of an experimental sequence we turn on the Doppler cooling laser and 
make a projective measurement of the ion spin state through state-dependent resonance 
fluorescence. With the Doppler cooling laser on, an ion in the |t) state scatters ~ 10 7 
photons/s, while an ion in \\) is dark. For the spin precession measurements summarized 
here we performed a global spin-state detection. Specifically we detected, with f/5 light 
collection optics and a photomultiplier tube, the resonance fluorescence from all the ions 



in a direction perpendicular to the magnetic field. 

We engineer an effective Ising interaction with a spin-dependent optical dipole force 
(ODF) obtained from off-resonant laser light. The ODF is directed either towards or 
away from regions of high laser intensity, depending on the frequency of the laser light 
relative to an atomic transition frequency. An optical dipole force is the analog of the 
well known ponderomotive force in plasma physics, here applied to the bound valence 
electron of the 9 Be + ion. We employ a spin-dependent ODF generated from the interfer- 
ence of two off-resonant laser beams at the ion cloud position, as depicted schematically 
in Fig. 3(c). If the frequencies of the two laser beams are the same, a one-dimensional 
optical lattice (i.e., ID intensity standing wave) is generated at the intersection of the two 
off -resonant laser beams. With a frequency difference /j,r between the two laser beams, 
the intensity standing wave becomes a "moving wave" with the intensity maxima pass- 
ing through the single ion plane at a frequency fiR. This generates an oscillating optical 
dipole force oc cos(jU^) . By adjusting the frequency and polarizations of the laser beams, 
we generate forces that are equal and opposite on the ||) , ||) spin states [11], resulting 
in the ODF interaction, 
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Hodf = - £ ^ocos (n R t)zj ■ &) ■ (1) 
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Here N is the total number of spins in the array, Zj is the axial position operator for ion 
j, and dj is the z-component of the Pauli spin matrix for ion j. A typical optical dipole 

force Fq used in our work is ~ 2 • 10~ 23 N. 

As discussed in the previous section, Hodf can generate an effective spin-spin in- 
teraction. Quantum mechanically, one solves for the evolution operator associated with 
the Hamiltonian Hodf - Because Hodf explicitly depends on time, this could involve an 
infinite sequence of commutators of Hodf with itself at different times. However, this 
sequence truncates after two terms. One term involves a product of Zj and &j, and de- 
scribes spin-motion entanglement generated by Hodf [20]. This term can be minimized 
by adiabatically turning on and off the interaction, as discussed in the simple examples 
in the previous section. The second term, which we label Hj, involves the product of 
different spin operators and describes the effective spin-spin interaction [21], 

Here Jjj is the strength of the effective Ising interaction between ion (or spin) i and ion 
(or spin) j. The // /s can be expressed in terms of the transverse (or drumhead) modes 

^COm^bmj of the array, where co m is the mode eigenfrequency and bj_ m (the jth element 
of b m ) denotes the relative amplitude of ion j, 
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FIGURE 4. Explicit calculation of the pairwise coupling coefficients J, j for N — 211 spins and plot as 
a function of spin-spin separation, du. For — co z < 1 kHz, Hodf principally excites COM motion 
in which all ions equally participate: the spin-spin interaction is spatially uniform. As the detuning 
is increased, modes of higher spatial frequency participate in the interaction and j develops a finite 
interaction length. 



Equation 3, along with a calculation of the drumhead modes, can be used to determine 
the effective spin-spin interactions. Figure 4 shows the results of such a calculation for 
different positive detunings of the optical dipole force beatnote /Ir from the center-of- 
mass (COM) mode 0) z . Note that the dispersion of the drumhead modes is anomalous; 
the longest- wavelength mode (co z ) is the highest-frequency drumhead mode. Figure 4 
shows that for ji R > co z , the pairwise coupling coefficients are positive (/, j > 0) , result- 
ing in an antiferromagnetic interaction, in agreement with the simple arguments of the 
previous section. As the detuning /i^ — co z is increased, the overall interaction strength 
weakens, and the range of the interaction gets shorter. In general the range of the interac- 
tion can be characterized by an approximate power law 7 ;J ~ // |r; — rj\ , where a can 
be tuned between and 3. The spin-spin interactions readily generated with trapped ions 
are long range. These long-range spin-spin interactions are "unphysical" in the sense 
that the strongest spin- spin interactions in solid-state materials are due to the exchange 
interaction, and therefore nearest neighbor. However, long-range interactions provide a 
good mechanism for generating long-range quantum entanglement, which appears to 
be a property of states with nontrivial topological/quantum order [22]. Therefore the 
long-range interactions readily obtained with trapped ions provide an interesting con- 
trast with naturally occurring solid-state systems, and an interesting place to search for 
new emergent behavior. 

The assumptions utilized in deriving Eq. 3 are substantial. For example, we assume 
that each ion satisfies a Lamb-Dicke confinement criterion [11] based on the "wave- 
length" of the optical dipole force, and that the optical dipole force wavefronts are pre- 
cisely aligned with the ion planar array, to name just a few. It is important to benchmark 
the experimental system to assure that it behaves like Eqs. 2 and 3. We do this by com- 
paring experimental data with the mean-field prediction that the influence of Hj on spin j 
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FIGURE 5. Measurements of /, normalized to the measured laser intensity, as a function the ODF 
beatnote detuning from the COM mode. The solid line (red) is the prediction of mean-field theory that 
accounts for couplings to all N transverse modes; there are no free parameters. The line's breadth reflects 
experimental uncertainty in the angle Or — 4.8 ± 0.25°. The mean-field prediction for the average value 
of the power-law exponent, a, is drawn in green (right axis, linear scale). 

can be modeled as a magnetic field Bj = jj Y,fy=jJi,j{&f) m me z direction due to the re- 
maining N —\ spins. The details are discussed in Ref. [1 1] and supplementary material. 
This interaction-induced magnetic field gives rise to spin precession in excess of normal 
Larmor precession. We prepare all the spins in the same initial state, rotated by an angle 
9 from the z axis (i.e., prepare spin \\) and then rotate by 0), and measure the excess 
spin precession as a function of 0. Through global fluorescence measurements we mea- 
sure the excess precession averaged over all ions, from which we are able to extract the 
average pairwise coupling /= -^j Ylj=iHfy=jJi,j- Figure 5 summarizes measurements 
of / as a function of the beatnote detuning /Ir — (O m . With no fitted parameters these 
measurements agree well with the predictions of Eq. 3, including contributions from 
all modes. Also shown in Fig. 5 is a theoretical calculation of the power-law exponent, 
which describes the range of the interaction. At the mean-field level, our system is well 
described by Eqs. 2 and 3. 

Future efforts include benchmarking quantum effects such as depolarization [23] and 
squeezing/anti-squeezing due to the quantum many-body interaction. An intractable 
simulation requires adding a non-commuting term to Hj, for example a transverse 
magnetic field, which we can do with the 124 GHz microwaves. Several technical 
improvements likely need to be made before this step. For example, the coherence of Hj 
is currently limited by spontaneous emission from the off -resonance laser beams. This 
can be greatly reduced by increasing the angle Or between the two ODF laser beams. 
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